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STOPS IS 

Shallow footings are the oommonly occurring typg of 
foundations for most Civil Engineering structures., The 
general practice of footing design is first to estimate the 
plan dimensions- providing adequate safety against bearing 
capacity failure only, and then checking for settlement; if 
the settlement is with in the permissible value, detailed 
structurial design are carried out on the basis of the obtained 
footing dimensions. Detailed investigations on the choice 
of proper shape and dimensions on the economy are seldom 
made. 

As such, in the present investigation an attempt has 
been made to provide an easy and a direct e^proach for the 
optimal limit state design of an isolated column footing, 
with due attention shown towards the settlement and other 
essential aspects for a safe design. Plan dimensioning of 
a group of footings is generally done by the conventional 
method based either on equal bearing pressure or on an equal 
settlement, A method of plan dimensionjng based on allowable 
differential settlement can be successfully applied to achieve 
a considerable saving in the plan area of the footing. Hence 
a generalised approach for the plan dimensioning of a group of 
footings has also been carried out. 
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The efficienoy •f numerical methods is problem oriented; 
so when they are applied to a new class of problems, there is 
a need for critical eyaluation regarding their applicability. 
Such studies have also been undertaken. Srom the investigations 
undertaken it has been concluded that penalty function 
technique along with Powell's algorithm for unconstrained 
minimization and quadratic fit for linear minimization are 
quite efficient for the problems that have been solved, 

Purther it has been observed that the developed 
general methodology to obtain the optim?il solution leads to 
a quite substantial amount of net saving over the conventional 
manual design. The net saving is about 10 to 40 percent 
depending on the type of problem. 




OHAPTBR - I 


INTRODUG'DIOU 

1,1 GSMiL ; 

Most of the time inadequate attention is paid to the 
economic aspects in the design of an engineering system. Biis 
results in the wastage of considerable amount of scarce 
resources which could have been saved and utilised in a better 
manner. As such, it is necessary to consider several alteranative 
designs in order to arrive at an optimal solution. Dais 
requires the use of computers to enable the designer to speed 
up the tedious and monotonous repeated computational work. The 
computer aided design of. an engineering system requires the 
implementation of sophisticated algorithms for the , analysis to 
minimize the wastage of computer time. 

Most of the real life problems of design and analysis 
are nonlinear in nature. The mathematical problems that arise 
in their study stretch the limits of conventional analyses and 
requires new methods for their successful treatment. In the 
use of analytic or numerical techniques for solving complex 
problems, the focus for discussion falls not only on the 
veorious techniques available for the analyses but also on the 
art of how such mathematical procedures are applied. Mo one 
procedure or series of procedures will be the panacea that solves 



Q-ll problems to the last detail. While optimization theory 
is well established for well defined systems and models, its 
efficiency is problem oriented. So in the study reported here, 
the usefulness of such techniques in the settlement controlled 
optimal design of isolated column footing and plan dimensioning 
of a group of footings have been presented. Parametric studies 
have been undertaken and reported. 

Design of an isolated footing involves the proportioning 
and the structural design of the footing. The structural design 
of isolated footings has received sufficient attention; but the 
first part of the design i.e., proportioning of the footing 
in contrast, has not received the same attention. 

The general practice of footing design is first to 
estimate the plan dimensions providing adequate safety against 
bearing capacity failure only, and carrying out the structural 
design on the basis of obtained footing dimensions. Detailed 
investigations on the choice of proper shape and dimensions 
are seldom made. It is a well known fact that the basic 
criterion governing the design of a foundation is that the 
settlement must not exceed specified permissible value. OOhis 
value will vary from structure to structure. My design which 
does not take in to account the settlement aspect in the 
analysis is a suspect from the safety point of view. However, 
in normal practice the settlement aspects are indirectly taken 
care of by using an excessively high factor of safety against 
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bearing capacity failure, but it results in an uneconomic 
design* 

The reported study presents an easy and direct approach 
of arriving at an optimal cost design of an isolated colimin 
footing satisfying the various foundation and structural design 
re quir em ent s . 

The bearing capacity and settlement characteristics 
of a group of footings is different from that of a single 
isolated footing because of the interference amongst the 
adjacent footings. The phenomenon of interference is of 
great importance in closely built in areas where there will be 
interference of stress distribution in the foundation soils. 

It has been investigated that one of the various factors 
affecting the bearing capacity and settlement of a group of 
foundations is the spacing between the foundations. 'The 
interference of footings on sand has been observed to cause an 
increase in bearing capacity fand decrease in settlement with 
the reduction in spacing, (Singh, et.al. 1973). For a safe and 
economic design, it is essential that the effect of interference 
on the stresses and displacements in the soil is not ignored. 

As such, a generalised approach for the optimal plan 
dimensioning of a group of footings has also been included in 
the reported investigation.. 



It is a well known fact that the efficiency of 
numerical procedures are to a great extent problem oriented. 

As such, even when an established good numerical scheme is 
applied to new class of problems, it needs critical 
evaluation. In the reported investigation, sequential 
unconstrained minimization technique (SMi) is used to arrive 
at the optimum solution. Interior penalty techniques along 
with the Powell's method of unconstrained minimization and 
quadratic fit for linear minimization has been used. Sffect 
of starting point, penalty parameter etc. on the efficiency 
of the numerical algorithm when applied to the enunciated 
problems have been studied and reported. 

1 . 2 -BRIli!? amiTUES RSYIHW; 

Very few literature is available on the optimal design 
of footings, most of the literature pertains to the structural 
aspect only with complete disregard to the settlement considers^- 
tions. However, there are plenty of literatures available 
on the various aspects of determination of stress distribution 
in the soil medim, bearing capacity and settlement analysis 
of foundations. Only those papers which are directly relevant 
to the present study are referred, 

1 . 2.1 Optimal Design of Footings; 

Optimum design of isolated column footing using 
sequential linear frogranming was candied out by (Bhavikattl, 
et.al. 1979), In this paper, the progress of optimization was 
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discussed and several parametric studies were carried out. 

It was concluded that a designer could achieve 8-10 percent 
economy by keeping the effective depth about 15 percent less 
than that obtained in the balanced section design. This 
paper does not consider the settlement aspects in the 
design. 

Subbarao, et.al.(l972) have given a method for an 
economic plan-dimensioning of footings subjected to uniaxial 
moment. However, the method does not take in to care of the 
tilt in the footing. White, 1956 has given a graphical method 
to arrive at the length and breadth of a rectangular footing 
subject to uni-axial moment. But, both Subbarao and White 
have not taken in to account of the effect of depth of 
embedment on the required plan-domensions of a footing. 

1 . 2. 2 Bearing Capacity: 

There are plenty of papers available on this aspect of 
foundation engineering. However for the sake of space and 
brevity, only those papers which are useful to the present 
work are presented. 

Meyerhof ( 1953) showed that the footings subjected to 
eccentric loads could be safely handled by introducing a 
fictitious effective width, 3' = 3-2e of the footing instead 
of the actual width, ’e' being the eccentricity. 



Significant contributions to the hearing capacity- 
analyses of layered soils have been made by Button (1953) > 
Sucklje (1954), Sivareddy and Srinivasan (1967), James et.al. 
(1969). Brown and Meyerhof (1969). 

The studies of Sivareddy and Srinivasan (1967), 

James et.al. (1969), both representing the extentions of 
Button's work, demonstrate the fact that the case of aniso- 
tropic layered clays and the case of single clay layer with 
continuously variable strength can be handled with sufficient 
accuracy by introducing an average strength for the layer or 

f 

soil zone in question. 

1.2.3 Stresses and Bis-placements : 

Yqvj frequently, we come across the case of a 
reotangular area subjected to a load of linearly increasing 
intensity, 2his case was solved in a comprehei^ive way by 
(xiroud (I 969 , 1970a, 1973), the results of his studies are 
set in tables and therefore convinient for practical use. 

But Giroud's procedure can not be employed when the 
load is not dis-feributed linearly tuid the loading area is not 
rectangular. In such cases one c^m proceed using Saint Tentmt' 
principle in the following way. 

The iictual loaded area is di-vided into a number of 
small, usually square areas, and the load applied to them is 
replaced by concentrated loads . The resulting stress is then 
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the sum '>£ the stresses produced by these concentrated loads. 

To attain an accuracy more than 95 percent at a given depth, 
the dimensions of the areas (squares in the case considered) 
must be smaller than 1/3 of the depth (Peda, 1978). 

Stresses and displacements in the dase of non- 
homogeneous soils were studied by G-ibson (1967, 1973). He 
presented an equation for the corner settlement of a 
rectangular area subjected to a uniform pressure intensity. 

Thenn - De Bjarro^s (1966) handled a three layered soil 
by determining a single composite modulus for the entire soil. 
Their method can also be extended to any nimaber of layers. 

"I • 5 MOTIVATIOIT AM) S OOP 3 0? TEJS WORK; 

It is a well known fact that settlement and tilt governs 
the design of a footing subjected to eccentric loads, especially 
in sands. Unfortunately, no work on the optimal design of 
eccentrically loaded footing has been done, which highlights 
the aspect of settlement on its design. 'The present work 
attempts to provide an easy and a direct approach for the 
optimal limit state design of an isolated column footing, 
with due attention shown towrards the settlement and other 
essential aspects for a safe design (Chapter 3). Tbe idea 
of an engineer should not only be on a safe design, but 
also on an economic design which fulfills the requirements 
of an ideal design. In this work, an attempt has been made 



to achieve economy in the design, satisfying all safety 
aspects. 

Plan dimensioning of a group of footings is generally 
done by the conventional method based either on an equal 
allowable bearing pressure or on equal settlement. 'Ihe 
conventional method of design leads to a lot of wastage of 
plan area of the footing. A method of plan dimensioning 
based on allowable differential settlement can be successfully 
applied to achieve a considerable saving in the plan area 
of the footing. Oomputer aided design can alone be a 
resort to this kind of problems involving several trial 
designs. 



CEAPTSR - II 


BASIC PRJ-EGIPISS Ai© BESIG-i'T PROOBDURBS ADOP’IHD IN TEh] PROPOSED S'lUDY 


2. 1 CENBRAI; 

The general principles used in investigating the 
problems as discussed earlier are presented in this chapter. 

The relevant basic concepts and approach of estimating the 
stresses at a point in the soil media, bearing capacity and 
settlement analysis of the fouridation and the limit state 
design are presented and discussed, 

2*2 LIMIT SPATE DESIG-N; 

It is generally agreed that the stress-strain behaviour 
of reinforced concrete (R. 0.) deviates from linearly elastic 
behaviour appreciably even at lower loads and radically at higher 
loads. Still, the 'straight line' behaviour of the stress- 
strain relationship or the linear elastic theory is being 
widely used in design becanse of its simplicity. Extensive 
research has conclusively established that lineiar elastic 
theory is inadequate for a comprehensive understanding of 
reinforced concrete structural elements under loads progressively 
increasing up to failure. Consequently, ultimate strength 
design (U. S.D. ) recognising non-x-linear stress-strain relation 
for concrete has been accepted as a means to determine the 
true ultimate strength of reinforced concrete sections. 
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All the ultimate strength design, studies on. reinforced 
concrete structures that have been carried out have essentially 
looked at the strength aspect only, pushing the serviceability 
requirement into the background. 'Ihe advent of ultimate 
strength design and high strength materials has resulted in 
more flexible structures for which the serviceability 
considerations are as important as the strength requirement, 
if not more. It is observed from practice that most common 
cases of damage in R. G. structures are excessive deflection, 
cracking etc., which affects the intended utility of the 
structure. The safety against these damages thus demand 
priority for consideration in the design. This has naturally 
led to the limit-state design philosophy. 

2.3 DERI ill TIOII 0? UMIT STA'TBS : 

A structure, or a part of a structure is rendered 
unfit for use when it reaches- a particular state, called a 
'limit state', in which it ceases to fulfil the function, 
or to satisfy the condition, for which it was designed. 

It is convenient to group limit states in to three 
major categories as follows; 

(i) Ultimate limit states; Those corresponding to 

maximim load carrying capacity and may be caused by the 
rupture of critical sections of the structure or excessive 
plastic deformation leading to collapse. 
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(ii) Serviceability limit states! Those corresponding 

to excessive deformation, wide cracking, excessive vibrations, 
undesirable damage etc. 

(iii ) Otner limit states; Those cor'responding to fatigue, 
ligh'tningjfire explosions etc. 

2.4 SJgBTY TAGTOR : 

The International Stadard Organisation (ISO), aiming 
at unification of different methods of structural calculations 
and ensuring the safety of structures, has recommended a semi- 
probabilistic limit state method. It is assumed that the, 
statistical distribution of loads and material strengths are 
available. ' Gharacteristic values’ of loads and material 
strengths are then selected with specified probabilities 
ag,ainst higher load or lower strength respectively. Ihe 
practice of using a global factor of safety has been replaced 
by a probability-based partial factor of safety format. The 
design loads are then computed by multiplying the characteris- 
tic loads by specified safety factors which vary with the 
degree of seriousness of the particular limit state being 
reached, probabilities of two or more loads occuring together 
and the reliability of the structural theories being used. 
Similarly, the design stresses are obtained by dividing the 
characteristic strengths by another partial safety factor which 
accounts for the difference between the strengths of controlled 
specimen and the material in the real structure, and any 
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other possible but unpredictable reduction in strengths. 

With these modified loads and material strengths, safety of 
the structure is investigated for the relevant limit states. 

2 . 5 DBSIQ-Ni 0? SHORT QODIPRESSIO'N MM3SRS U>r.DMl .3IAXIAI. 

The limit state design of short members subjected 
to biaxial bending is carried out by using the standard design 
procedure given by Mallick and G-upta (1980). However, it is 
necessary to discuss about uniaxial bending before proceeding 
for biaxial bending. 

Referring to Rig. 2.1, 


T 

= fst, Ast, 

G_^ = fsc. Asc 

(2.1) 

^c 

= 0.362 .fck 

. b . XU 

(2.2) 

P 

= ^c ^ <^sc - 

T 

(2.3) 


The moment capacity, MUXI = P.e, is given by: 

MUXI ^ P.e = 0^(0. 5 BG - 0.416 xu) + (Ogc+T) a 

(2.4) 

and a = (DC - 2 x cover) / 2 (2.5) 

where, 

T = force in tension steel 

Gq = compreshive force in concrete 

Ggg = force in compression steel 

fck = characteristic strength of concrete 

DO — depth of the column , 

P = axial load on the column 



13 


XU = depth of neutral axis 

b = breadth of column 

Asc = area of compression reinforcement 

Ast = area of tension reinforcement. 

The design of membeis subijected to biaxial bending is given 
by the equation, (I.S. 456-1978) as 
MUX a MUY a 

(« ) + ( ) < 1.0 , ( 2 . 6 ) 

MUX1 MUY1 

where MUX and MUY are the moment capacities about the x and 
y axes due to the design loads. MUX1., MUY1 are the maximum 
uni-axial moment capacities for an axial load Pu, bending 
about X and y axes respectively. 

Value of a depends on Pu/Puz 

where Puz = 0,45 fck . Ac + 0.75 fy. -^c (2.7) 

fy = Characteristic strength of compression 
reinforcement, 

Asc and .Ast are the areas of longitudinal compression 
and tension reinforcements respectively, 

= Area of concrete in the column. 

lastly, fsc and fst <ore the stresses in compression and 
tension steel respectively. 

The value of a is as follows: 

a = 1.0 for Pu/Puz <0,2 
and a = 1.0 - 2.0 for 0.2 < Pu/Puz < 0.8, 


H 


■Ihe critical section of rhe column is the one which 
is in contact with the footing. 

'Ihe column in the present study is short one and 
hence no end condition correction is necessary. 

2.6 SESI &! ! 0? gOQfllia; 

The following are the notations used in the design of 
the footing, (?ig. 2.2). 


W1 = inclined load on the column 
W1 can be resolved in to horizontal and vertical components. 
H1 = horizontal load along x axis 

H2 = horizontal load along y axis 

HX = Bccentricity in the 'x' direction 

EY = Bccentricity in the '3/ ' direction 

Df = depth of embedment 

= partial safety factor, equal to 1.5 
d = effective depth of the footing at the centre 
I) = overall depth of the footing at the centre. 

Ref erri ng 1 0 ( Bi g . 2 . 2 ) , 

and M^y- are the moments about x and y axis 
respectively, which also includes the moments caused by 
horizontal load’s H2 and Hi on them. 



i.e. 


WUX + H2 (Df - E) 
MUY + HI (Ef - D) 


( 2 . 8 ) 

(2.9) 



and 


The other notation are: 

B = breadth of the footing 
q-] = maximum soil pressure along x axis 

^2 minimum soil pressure along x axis 

= maximum soil pressure along y axis 

= minimum soil press’ure along Y axis, 

jd = lever arm. 


Assuming the self weights of the column and the 
footing as 5 percent and 10 percent of the axial load coming 
on the column, the maximum ;and minimum pross-ure on the soil 
are given by 




ax 


1.155 W 6.0 X Mux 

"f '■ 

IS 1^2 

6 • 0 X 


( 2 . 10 ) 


1.155 W 


ani 


'^in 


IlB 


,, 6.0 X Mujj. 

2 

IB 


6.0 X I'l, 


ay 


( 2 . 11 ) 


31'" 


respectively. 

Soil pressures, q., and q. 2 , q^ and are given by, 
(Big. 2.2). 


1.155^'' 6.0 X M, 


l-l = 

0.2 ~ 


. IB 

1.155W 


+ 


ux 


Ib2 
6.0 X 


(2.12a) 

(2.12b) 


IS 







1. 15515 6.0 X HL^ 


IS 


3L 


1.155 w 6.0 X Ivl^y- 

IS 3L^ 

iow the eccentricities, EX and ]i!Y are given by, 


(2.12c) 
(2. 12d) 


BX = 


and EY 
further, 

Pi 

P2 


1.05 W 
1 .05 w 


(2.13a) 

(2.13b) 


1.05 w 6.0 X ^ 


A 

1.05 W 


IS^ 




^ ^ ^ux _ £J!h2L. 


(2.14) 

(2.15) 


A 13*^ 

where, p-j and P2 pressure causing the bending moment on 

the footing slab and A is the area of the footing. 

Maximum bending moment at GO is given by, (7ig.2.2) 


3-b 


_B-b p 


where, P2+25: 
and z 


Ml =1^^ ^ ^2.d6)_ 

P2(3-b) + P-, (3+b) 


(2.17a) 


23 

= p.| - (p2+x) (2.17b) 

Ihe maximum bending moment about AA is given by, (Big. 2.2), 


1 • 1 

^2 “ ^ f ^ (p2+x) (S-b) ^ B + —— ( ZZ . (Ii-b) 3 ) 

u- 


(2.18) 
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7'(P1+P?) + b (p-t-Pp) 

where, Pp+x = 1 . f; (2.19) 

2L 

and ZZ = - (pp+x) (2.20) 

Further, the moment capacity of a trapezoidal section is 
given by 

MR = (( k-kp) kp ) d^ fck (2.21) 

where, b-| and bp are the widths of the section at rhe 
compression face and at the level of tension steel respectively. 


k and. kp are the coefficients whose value are equal 
to 0,145 and 0,03 respectively for concrete. 


The moment capacity MR in each direction is equated 
to the corresponding maximum bending moments M-] and Mp 
respectively in that direction, and the effective depth ’d' 
is determined in both the cases by the equation, 

jMaximum bending moment 

d =f :: (2,22) 

MR 

The higher value of 'd ' is selected as the effective 
depth 'd' , which is further checked so as to ensure that; it 
satisfies the requirements of one-way and two-way shear. 


lastly, 

and breadth (3) o 
AST2 = 


area of steel required along the length (Ii) 
f the footing are given by. 


1. 15 Mp 1.15 X 

and .AST3 = — 


(2.23) 


D-d-.fy 




respectively 



2.7 DSSI&jI 0EI3GSS J-OR S.OR’IY: 


The following design restrictions are checked for 
Scifety, as described below. 

2*7 .I S i ag 0 n al Ten si o n : 

'The nominal shear stress at the section 'aaaa' 

Rig. 2.2, is given by 

“ “d" 

0 

1.05 W . ^ .2s . . 

where, Y = '0 ^ CSS — (b-t* d) ) is the shear force at 

the section and b^ = 4 (b+d) is the peripheral distance of the 
critical shein: zone. 

The allowable shear stress CZa) ^^15 concrete given 

by» 

^^=0.25THk (2.25) 

Ror safety against she:ar failure, 

Lv 1 Ca- (2.26) 

Rotations used are the same as those described in the begining 
is chapter. 

' .C^ eck for O n e-Way Shear ; 

Ref er ’ing to Rig. (2.2), the shear force (s) at the 

cal section '00' is given by, 

3-b 

S = p.l (( — ) _ d) (2.27) 

, b = breadth of the column, 

d = effective depth of the footing, 
p = average soil pressure. 



-forainal shear stress (q^) at rhe critical section '00* is 
given by 

S 

^ = (2.28) 

b 'd ' 

where, b' = width of rhe footing at the critical section 
and d' = effective depth at the critical section. 

Hie safety against one-way shear, 

< % (2.29) 

where, is the allowable shear stress as given in 

I.S. 456-1978. 

Ihe ■‘■hicknesa of the footing provided at any point 
is also checked for adequacy. 

2.7.3 Check for Searing Stress (f^'dj:)* 

Searing stress calculated at the interface 

between the column and the footing. 

fhe IS 456-1978 has specified the permissible bearing 
stress depending on the characteristic strength of concrete, 
supporting area of the footing and loaded area of the column 
For safety rhe actual bearing stress should be less then 
the permissible value. 

©^r < 0.45 fck ^ A^7 A 2 

TA 1 /A 2 <2.0 


with 


(2.30) 

(2.31) 



wnere, 


- supporting area for bearing of a footing and 
-^2 = loaded area at the column base. 

2.7.4 f actor of Safety Against Sliding; 

lactor of safet 3 r against sliding (?.S.) along a 
direction is given by, 

^L^ Y, 

P.S. = — — (2.52) 

S H 

where, V = lotal weight of the isolated column footing 

including the backfill and the load on the column, 
■iL = coefficient of firction between concrete and 
the soil, 

and 2H = sum of the norizontal forces in that direction. 


In the present case, H1 and H2 are the only horizontal 
forces in tneir respective dircections, 

lor safety against sliding, I’.S. > 1.5. 

2.7.5 Check for Overturning (S.?) : 

It is the ratio of restoring moment to the moment causing 


overturning in that direction. 

.. . r. ~n _ Restoring Moment (9 

i.e. , S.R. = — — — — — — 

Ovc5u lurning Moment 


where, . 

Restoring moment about the length L of the 
footing . 

lurning moment about the length 1 of 
the footing 


= 7x3/2 

(2.33b) 

(2.33c) 


= HIxRf 



= 7 .X L/2 (2.33d) 


Simtl.3rly, restoring moment about 3 of 
the footing 

Turning moment about the breadth 3 

of the footing = H2 .x 33 (2.33c) 


Restoring moment is the moment of all the vertical 
forces acting on the footing about its ce^it^oid and turning 
moment is the moment caused due to the horizontal load plus 
the moment acting in that direction. 

2.3 G.giOrTY; 

The general aspects of bearing capacity are described. 

below: 


factors that determine the ultimate be;aring capacity 
of horizontal footing with a central load include (1) unit 
weight, shear strength and deformation characteristics of 
the soil, (2) the size shape , depth and roughness of the 
footing, ( 3 ) the water table conditions and initial stresses 
in the foundation soil. 


Hansen 


The general bearing capacity equation suggested by 
(I970)is given by. 




ulr 


Q ^^c ^c ^c 


3^ Sq, d^ i^ + 0.5^) BNS id 


^ V 
(2.34) 


where, G = unit cohesion 

■' , , [ 

3 = least lateral dimension of the footing 

.3 , 3 and 3.1 are the Ter zaghi 's bearing capacity fao-cors| 

C % Y ■ I 
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S^jSq^ and 3^ are the shape factors. 

Iq , iq^ -and are the lo?ad inclination factors, 
lastly, dg, dq and d,^ -are the depth facnors. 


iiet bearing press lare is given by = '^ult' ~V ^ % 

where lu.lt' “ gross bearing pressure (2.3'' ) 

Ihe allowable bearing pressure (ly^) is the net taejaring 
pressure divided by an appropriate safety factor (iF) . 

„ = (2.35) 

where, '0^ = depth of embedment, . 


2 . S . 1 Wa-cer Tabl e lilffec ts: 


Ihe presence of water table around a footing reduces 
the effective shear strength of a gr annular soil and hence, 
its bearing capacity is also reduced. For a fully submerged 
footing', ,the unit weights used with the and term..s of 
the bearing capacity equation is the buoyant unit weight. 

Since the buoyant, unit weight is about one half of the moist 
unit; weight, the bearing capacity of a submerged footing is | 
about one half of that of the footing well above the water 
table, if the water table is at a’ di stance 3 below the 
footing, it' is'' Assumed t:b have no effect. Tnen the %a ter 
table is at tne base of the footing, the buoyant unit weight 
is used only with the, 11 term. 

For water table positions intermediate between the 
g^imd-suffhce' of the footing or between the I 



tiase of the footing and a distance .3 below the footing, 
the respectiYe unit weight is adjusted by linear interpola- 
tion. Then method is available in any standard text book 
on foundation engineering, (Bowles, 1977). 


2 . 8.2 Ec c en tr i c lo ads : 

footings, must frequently resist an overturning moitrent, 
which results in an eccentrically applied load. lieyerhof 
( 1953 ) has proposed that the dimensions 3* and 1 ' of an 
equiv;alent conccentric footing be found and used to determine 
the allowable bearing pressure. 


The equations for 3' and L' are 
B' =3-2 e -g and 

jjf — Xi — 2. 


( 2 . 36 a) 

( 2 . 36 b) 


where 63 and e^ are the eccentricity in each respective 
direction. There total allowable load is qj^(3'L*). 

2.8.3 Bearing Capacity in Lay ered Soilsl 

Iwo oases ot bearing capacity analysis oi layered soils 
are considered. They are , 

Case 1. Soft layer overlying a stiff layer, 3ig. (2.3a) 

Oase 2. Stiff layer overlying a soft layer, 3ig. ( 2 .^a) 

The bearing capacity of the foundation in both th 
cases is given' by '.Brown and Meyerhof (1969), as 


(2.37) 



v'/here G1 = undrained shear strength of the upper layer ;md 
a modified be,n:ing capacity factor which depends on, 
relative shear strengths of the two layers, K =C2/01, the 
relative thickness of the upper layer, H/B, as well as the 
foundation shape. 

Qq = ultimate bearing capacity 
Vesic ( 1970 ) suggested the following expresvsion for the 


modified bearing capacity 5^^ for the situation of a soft 


clay layer overlying a stiff clay layer. 

(yj + P-OQk+I) (1+Kp ) + ^- 1 j 


% 


[K(K+1) N*+K^_i2 ((FqHp) 1 * + (K.lJ+ P-1)(]75+1) 

( 2 . 38 ) 

B1 

where p = 


* 


2(B+L) H 


B = breadth of footing 

1 = length of footing 

H = thickness of the upper layer. 


Bor the case of stiff clay layer overlying a soft cla^ 
l;y/-er, Brown and Meyerhof (1969) suggest the following 
expression for 


-7, 


‘M 


+ K 


C ^'^0 


^ ii'a 


(2.40) 


where Q is the correction factor. | 

A more general analysis has been given by (Yesic, l970) 
for a rectangular footing resting on an upper stronger layer 
having strength parameters Cl and 01 underlain by a lower lay^ 





of strength parameters 02, and jZS2. The expression for ultimai; 
bearing capacity (q^) is given by 

^ ^ ..... ^ -j 

'lo " v^o" + ^ [2(1+3/!.) K tan ^^(e/B)\ 

1 

- X O-i X Got 0., (2.41 ) ' 

I 

where K = (l _ Sin^ 0^) / (1 + Sin^ 0^) (2.42) 

3Jad q^’ = bearing capacity of the footing resting on the 

top of layer 2, to be evaluated with strength parameters C 2 an 
^2 and other characteristics of the second layer. 

2 . 8.4 Procedure for Oalculating the Average Value of 14, Simso 

(1977): 


(a) 

(b) 

(c) 


The average value of d is calculated as follows : 

The average value of V between the foundation level i 
and depth of 3B/4 below the foundation is taken and Is! 
multiplied by three giving 3.^71 . 1 

Average value of d between depths 33/4 and 33/2 is | 
tiaken and is multiplied by two giving 2d2. 1 

lastly, the average value of 1 between 33/2 and 2.3 

( 

is found, ( 33 ). j 



( 2 . 43 ) 


waere, d is the standiurd penetration number 



^ ® • 5 P rediction of Allowable Bearing; Pressure ; 

Peraaghi and Peck have given in the following 
e-xpresaion for allowable bearing capacity for 2.5 cm settlement 


of individual footing, differential settlement of 2 cm. 

lall =3.5 (.V3) Ra (2.44) 


wnere- 


3.5 (.V3) R,., 


~ allowable net increase in soil pressure over 
existing soil pressure for a settlement of. 

p 

2. 5 cm in t/ra . 

= Standcird penetration ntimber with applicable 
corrections. 

3^^ = lldth of footing (or least lateral dimension) 

in meters. 

= Depth factor = 1+ 

R(i < 2.0 (2.45) 

R^^ = Water reduction factor 

= 0.5 ( 1 + %/3o) (2.46) 

= Depth of water table below foundation 


(2.47) 



At 

% 


0, 

Ryi^ 0 . 5 

and 

.A.t 

% 

= 

% » 



Meyerhof (1956) has proposed slightly different equations, but 
they yield for the usual size of the footings almost the 
same results. 

Ihe aoove equation is considered conservative. ThG 
net pressure may be increased by about 50 percent (.'Bowles, 1 968) 




2* 9 ,gOjCT,TATIO.J O'? SJRESSlilS IK THS SOIL M.ASS : 

'^ery frequently, we come across a rectangular area 
subject to a load of linearly increasing intensity. ‘Ihis 
case was solved in a comprehensive way by iliroud (1969,1970a, 
1973); the results of his studies are set up in tables and 
therefore convinient for practical use. 

'Ihe other, sufficiently reliable alternative could 
be to divide the actual area in to a number of small, usually 
square areas, and the load applied to them is replaced by 
concentrated loads, lo attain an accuracy more than 95 percent 
at a given depth, the dimensions of the given area (squares, 
in the case being considered) must be smaller than 1/3 of 
the deptn ( Pavlor 1948; llorin 1959; Myslivec al, 1970, etc.). 



Q = vertical load (concentrated load) 
r = radial distance of the point 
2 = depth below the ground surface 

~ vertical stress in the soil. 
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2.10 OOi^U PATIO 111 O'T S IT : 


2.10.1 Tne Buisman - Debeer Method: 


This method, in which the settlement under a given 
pressure is estimated, depends on the correlations between the 
cone penetration resistance and the corapressibilit;/ of sand. 
Buisman proposed the following emperical equation to obtain 
the value of compressibility (O) from the value of the cone 


penetration resistance (q^) 


[ 2 . 50 ) 


where 


^ • b Iq 

crtr ‘ 

is the effective overburden pressure at the depth 


of measurement. 


Settlement can be obtained by means of the equation; 


H 


S = 


In 


( 


(To + A.O- 




(fi.5l ) 


0 crs 

where, S is the settlement of a layer of thickness 
H and A.x-is the increment of vertical stress at the centre 
of the layer. If 0 represents the compressibilitv of sand 
over an elemental layer of thickness ^ 3, then the settlement 
can be expressed as 


/ 

0 


H 1 


»' (in (“ 


+ 


) ) dz 


( 2 . 52 ) 


or approximately, 

H 

1-53 


3=2 ^ - 3-,: 


- ! 


' t , * — t- 

,v /-a \ 

A^lOg^Q 


(2.53) 
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■Based on a study of case reco37ds, Meyerhof (1965) , 
recommended tiiat the foundation pressure producing the 
allowable settlement by the 3uisman-Debeer 's method should be 
increased by 50 percent. This is approximately equal to 
using the following equation for the constant of compressibility 


C = 1.9 q^ (2.54) 

Schmertmanii's metiiod (1970) could not be used owing to the 
fact that the equation given by him holds good only when 
the area is subject to an uniform pressure of intensity q, 

2.10.2 .(^n solidation Settltanent ; 


In order to estimate consolidation settlement, the 
value of either the coefficient of volume change or the 
compression index is required. 

The settlement of a layer of thickness 'H' is given by 


H 


I 


Sc = / mv /\cr dz 
o 


(2.55) 

If mv and are assumed constant with depth, then, 

Sc = mv Air-H (2.56) 

3q _ a. (2.57) 


or 


^0 ~ H 
1 + e^ 


or, in the case of normally consolidated clay: 


Sc. = 


H G_: log (“rh-) 

^ ^ V r> 

4 + 


(2.58) 
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where, 


= initial Yoid ratio 
e-| = final void ratio 

= increase in soil preseore due to the load 


mv = 




..-3 


,TT 


Gompression index of the soil 
coefficient of volume change 
overburden pressu.re at any depth and 


+ f 


“S' 


(2.59) 


2 . 10.3 Settlement in the Oase of layered Soils : 


An approximate analysis can be undertaken by a simple 
extrapolation proposed by Steinbrenner ( 1934 ). Steinbrenner 
suggested that the displacements^,\'g_ of the point A on the 

J 

atnoface of limited depth S is computed from the Boussinesq's 
expressions for the displ^aoements of the points A and 3 
and%^ 3 ') for an infinite depth of soil (3ig. 2.3b). Since 
the deflection of the point 3 is zero,/h.^ A is given by 




(2.60) 


where, Zhij'A' and/g^-’g’ are the Boussinesq’s displacements in 
infinite depth of soil. 


A mul'cilayered system as in Big, 2. 5b is considered 
and the displacements ^ , a 52 '< 1 j 3 ’ an infinite 

depth of soil of modulus 3 (i.e. , as if the whole half-space 
had the modulus of the upper layer) is calculated, liie 
settlement in each layer is given by 

= A5t •" 6J2 
^)2 2 " ■ 


(2.61) 

(2.62) 
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=^53 -^^4 


( 2 . 63 : 


3ut all the settlements are based upon the use of a uniform 
value of 3 = El .'fo accomodate variations in E with depth, 
we snould therefore multiply the nth. layer of El/^* Hence 


the total settlene.nt is given by 


-j— etc. 


(2.64) 


Bor the case where 'B' increases linearly with depth, the soil 


can conviniently be divided into finite layers, of equal 


thickness and ascribed an average B for each, etc. 


'Phe rate of increase of ?3 with depth is expressed in terms of 


B_ (modulus at the surface). Phus, 


( 1 + s - 


(2.65) 


I';iitial settlement for a soil can be computed by using undrainedi 


modulus 'B^' and Poisson's ratio. ' aX,' in ihe ^Soussenesq 's 
expression for displacement (w) gi■^^en b^r. 


Q (1+U.) j Z" 


2(1 -.A) 


w = 


R J 


(2.66) 


wnere ,Q 


concentrated load 


Z = depth at which the displacement is required 
R = radial distance of the point under consideration.. 
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2.10.4 .Al lowable Settl emen;^ 

Settlement is important, even iihough no rapture is 
imminent, for three reasons! appeartance of the structure; 
utility of the structure; and damage to the structure. 

If the maximum settlemejit, 

denotes the minimum settlement, the differential settlement 
is the larger settlement minus the vsmaller. Differential 
settlement is also chfuracterised by angular distortion b/l, . i 
which is the differential settlement between two points divided: 
by the horizontal distance between them. 

Dambe and Whitman ( 1969) have given that the allowable : 
differential settlement for an 1. 0. 0. building fr^ame as I 

0.0025-0.0041 > where is rhe distance between the adjacent | 
columns that settle by different amounts, or between any two [ 

' E 

points that settle differently. [ 

!; 

I 

Mac Donald and Skempton (1955), made a study of 98 ! 

i 

i 

buildings, most of rhem being R.0.0. structures. fney concluded 
that a differential settlement of about 3.75 cm on clays and i 
2.5 cm on sands between two isolated footings does not affect 
the satisf actor;^'' performance of the structure. However, while i 

i 

considering the interference between adjacent footings, a j 

differential settlement of only 2 cm has been allowed between 
any two foundations. i 
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2* '0*5 j^lP-Wable ^Ib tal Se ttlement; 

L;arabe and ^'■/lii'tenan,, (1969) have siven the following 
allowance for the total settlement: 

Settlement ;< 40 mm for footings on s.ands 

and Settlement < 65 mm for footings on clays. 

2.10.6 Influence ?aotor for footings at a Depth : 

fox ( 1943 ), proposed a correction to apply to the 
surface deflection which depends both on the depth of embedment 
and poisson's ratio 

fherefore, the corrected settlement S^ becomes, 

Sq = S. 53 (2.67) 

wnere, 1 -^ is obtained by programming fox-equations to obtain 
the settlement ratio vs (Df/3) and for selected (V-'B) and 
poisson's ratios 

S = surface settlement. 

Expression for 1 -^ can not be presented here as the expression 
is too lengtny to present. Expression is given by fox in 
his paper, ( 194S) • 

2- 1 1 EASIO GOd a EPTS 0? O P TIMIZAIIOE f:Ji<!0RY: 

2.11.1 &ejieral:^ 

Optimisation is the act of obtaining the best result 
under given circumstances and satisfying all limitations .and 
restrictions placed on it. In the analysis or desi^gn of any 
engineering system, the ultimate goal is to minimise the 
effort required or maximize the desired benefit. Since the 
effort required or the benefit desired in any situation can 
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be expressed as a function of certain decision variables, 
optimization can be viewed as the process of finding the 
conditions that gives the maximuiti or minimum value of the 
f'unction. dgain, since the maximum of a function can be 
found by seeking the minimum of the negative of the same 
function, optimization am be taken to m.ean minimization, 
without any loss of generality, 

2.1 1.2 jle rniinol ogy ; 

Before coming to the discussion of various techtiiques 
of optimization, it is important to define certain basic 
terminologies which are generally used in optimization theory, 

'Ihe numerical Quantities for which the values are to 
be chosen in producing a decision or a solution are called 
'decision variables’, A vector B containing all the decision 
variables in a particular problem, will often be referred 
to as a 'decision vector’. It should be emphasized at this 
point that 'a decision’ is simply a set of values of the 
decision variables, even if it is absurd or inadeguate in terms 
of functional requirements . If a decision meets all the 
requirements placed on it, will be called as an 'acceptable 
decision’. When a decision , fails to meet any of the require- 
ments, it will be called as an 'unacceptable decision’. 

Ihe restrictions that must be satisfied in. ox’der to 
produce an acceptable decision or solution are collectively 
called as 'constraints’. A coz^traint -that derives from the 
performance or behavio:? requirements is called a 'behaviour 
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constraint ' , 

Anv decision problem which uses optimization technique 
will ha^re it's objective to obtain the best decision. The 
function with respect to which the decision is optimized, is 
called the 'objective function', and is designated as J or 

fS). 

J'ow a standard form of problem statement is considered. 
Inis is presented in a gener;alized form and the actual physical 
situation need to be taken' into account while formulating parti- 
cul.ar problems. After making the appropriate engineering 
judgements and defining all the necessary functions and limita- 
tions, the optimization problem is stated as follows; 

Minimize f (D) ,3^ = (d^ d^ 6.^) ~ 

Subject to g^(3) _< 0 j = 1,2, m 

1^(3) = 0 j =1,2, k 

where, g. and 1. represent inequality and equality constr^3ints 

V 

on the decision variables, respectively. 

An optimization problem in which there is no restric- 
tion on the choice of decision vector, is called an uncon- 
strained optimization problem; otherwise, it is called a 
constrained optimization problem. Most of the practical 
application problems are ccmstraimed; but some of the most 

powerful and convenient methods of solving constrained 

0' 

problems involve the conversion of the problem to one of 
unconstrained minimization, iind this technique has been used 


in this thesis. 
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** ' • 3 .^.e . JJncp ns tr ained Minimi zation 'Je chni que : 

Ihe Tnultidimensional unconstrained minimization 
technique that will be used in this thesis is, Powell's (1964) ' 

conjugate direction method. Detailed descriptions of this 
method is obtainable in any standard text on optimization 
(e.g. , Pox 1971 , Rao 1978). In this method, at each point 
in the function space under consideration, there is a preferred 
direction along which the values of the decision variables are 
changed systematically by a well defined scheme, which leads 
to the minimum. 'Ihus starting from an inixial guess one 
arrives at the minimum through a sequence of improved 
approximations, each derived from the previous approximation. 

'The various raetnods follow the same principle stated above I 

but differ only in selecting the direction along which they 
se.arch for the minimum. Powell’s method is a nongrad.ient method. 

Powell’s method has the property of quadratic convergence. ^ 

I 

Without going into mathematical det;h.ls, a brief description | 
of the algorithm is presented as follows. 7 

I 

I 

2 • 1*4 Powe l l’s Conj ugate Direction Method: 

This method is an extension of the idea of pattern 

move, for an intutive understanding, Pox (1971) described 

' 

the method as follows: ’G-iven that the function has been 

, 

minimized once in each of the coordinate directions and then 

[ 

in the associated pattern direction, discard one of the 

■ ; i 

' . " I 

coordinate directions in favor of the pattern direotion for 
inclusion in the next m minimizations, since this is likely 
to be a better direction than the discarded coordinate direction. 
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■After the next circle of minimizations, generate a new pattern 
direction and again replace one of the coordinate directions', 
fhe flow di-agram shown in ?ig. (2,4-) expl.ains the algorithm. 
'.The basic method has a tendency to choose nearly dependent 
direction in .ill conditioned problems and the method may fail 
to Gonrerge to the actual minimum. One simple remedy is to 
reset the directions to the original coordinate vectors 
periodically and/or whenever there is some indication that the 
directions are no longer productive. Powell has recommended 
a more effective procedure to avoid this difficulty(.?ox, 1 971 ) . 

'Hhe search is terminated when the relative change in 
the function value .'and in all the decision variables, between 
two consecutive cycles of minimization, is less than the 
desired accuracy. 

2*11.5 Quadratic Interpolation; 

In this tecnniiue the function P(of') is approximated 
by a function H (a) which has an easily determinable minimum 
point. H(a) is expressed as, 

2 

H(oc) = a + ba + coc 

the minimum of which occurs where 
dli 

— — = b + 2ct^ = 0 
da 

b 

or a 

The constants a, b and ,c for the approximating quadratic can 
be determined by sampling the function at three different a 
values, a.| , and a^, and solving the equations,^^^^ ^^ 
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5'^ = a+ba^+ca^ 


and 


-2 


^3 


a + ttttg + ca^ 
a + boc^ + qo,^ 


where ?■] denotes i'(a^), and so on. Solving the above three 
equations, the values of a, b, c and are obtained ass, 

c = ( j?.! ) {p. ^—o. j?.j ) (a^— a,| ) / 

(aj-a^) (a^-a^) (a^-a^) 

b = (;?2 - ?^ ) / (a 2 -«-,) - c (a 2 ^^ 1 ) 

2 

a = - ba I - coc-j and a.* = _ b/2c 

iTow cnoosing 0, t and 2t for a-], 0^2 ^^*3. a^, one obtains, 

a = 

b = (45'2-3 ?i-:? 3) / 2t 

2 


and 


c = (3'3+:?^-25'2) / 2t' 
a-^= (4V2-3?i-5'3) t/(41V25'3-25’-| ) 


?or o£.’i<i‘ to correspond to a minimum of H(a) , it must satisfy 


d^H 


da' 


* 


a=a 


> 0 


Since H is a qudratic, this requires 
c > 0 


or 


) > 2? 


fhe point a* is considered to be a sufficiently good 
approximation to the minimum of F(a) if 
iHCa*)-;?^^*) j 
l?(a*) I 

!Bb.e flow diagram shown in Fig. 2.5 explains uhe algorithm. 
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^ 6 function IJetliod : 

Ihe basic ob;iect of Penalty function Method is to convert; 
the original cons trained problem into one of unconstrained 
minimisation by blending the constraints into a composite 
function and making it possible ro ignore them at the 
minimization stage. In this method, tne numerical solutions 
are sought by solving a sequence of unconstrained minimization 
problems. 

Ihe advantage of Penalty function method lies in the 
fact that the powerful, well studied and reliable algorithms for 
the unconstrained minimization of arbitrary functions can be 
used in this met nod. 

Ihe penalty function formulations for inequality 
constrained problems can be divided into two categories; 
interior and exterior. In the interior formulation, the 
unconstrained minima all lie in the feasible region and 
converge to the solution as a special parfsmeter is varied. 

In the exterior formulation, they lie in the infeasible region 
an converge to the solution from outside. Ihe advant<ige of 
interior penalty function method is that, given an initial 
acceptable decision, i'c produces an improving sequence of 
acceptable decisions. Moreover, we approach the constraints 
in such a way tnat they become critical only near the- end,. 

This is a desirable feature in engineering decision process 
because, instead of taking the optimum decision, one can choose 
a suboptimal but less critical decision if reqijired. After 
studying the relative advantages and disadvantages of the 
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available me'bhods, the interior penalty function method was 


selected for solving the problems involved in this thesis 


In the interior penalty function method a new function 
(jZ5-f unction) is constructed by augmenting a penalty terra to 
the objective function. 'Ihe j2^-function is defined as 


0 (5,r^^) 


( 0 ) 




m 1 
D=1 ) 


where ? is to be minimized over all D satisfying 


(3) < 0 , j =1,2, m . 

Ihe penal t]?’ p.iraraeter 'r' is made successively smaller in 
order to obtain the constrained mi ni mum of 1 . 


'Ihe flow diagram shEran in ?ig. 2.6 explains the 
algorithm. 


The interior penalty function method requires a 
feasible starting point for the search towards the optimum 
point. In cases where finding a feasible starting point is 
difficult, the penalt^r function method itself can be used 
for finding it (Rao, 1973). But the procedure involved is 
time consuming. 


A promising alternative approach is to use the extended 
penal 'cy function proposed by Kavlie (Basudhar, 1 976). In this 
method the infeasible starting points are acceptable to the 
minimization algorithm. In this approach the 0-function is 
defined as _ 1 

p(p) _ r ^ =r- ; g-i (2) < S 

J *3=1 Sj(D) 3 


0(D, ric) = / 


L^CB) - r^j. 2 , (23«g^(D))/3‘"; g^.(3) > 3 


D=t 



4-1 

where m is the total number of satisfied constraints and 1 
is the total number of violated constraints, 

’ = - r^/6^ 

6^ = a constant that defines the transition between 
the two types of penalty terms. 

fhe basic principles of the optimization techniques 
used in this thesis has been covered ver 3 ^ briefly, fhe detailed 
description ^and explanation of these methods are available in 
all standard texts (Fox 1971, Rao1973). 

Ihe basic principles described in this chapter will 
be referred to in the subsequent chapters. 

2.1 1.7 Ci onvergence Griterion; 

j3ae convergence is assumed to be achieved when the 
change in the values of ob.-jective function and the decision 
vector, between two consecutive cycles, is less than the 
desired accuracy, i.e., 

< s. (2.69) 

F„ "" 

9. 

and 1 Fq - 1 ^ (2.70) 

Ilhe same convergence criterion is used all throdgb-, wherever 
checking of covergence is involved. 
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CHIPrSR - III 


COiIROIJj'i;D OP IM il OOSI OlilSIGr.!;^ 

OP AI iSOIAriJP OOLTR.L'T POOIIIO 

3.1 M ^BR AL; 

In this chapter, a general approach for the settlement 
controlled optimal cost design of an isolated column footing 
has been presented. Ihe structural and foundation design 
aspects have already been discussed and presented in 
Chapter 2. As such, herein, only the objective .function, 
design variables and the design constraints are identified, 
whenever found necessary references have been made to the 
appropriate sections of Chapter 2 for further details. 

5, 2 -rOIAL QUOIUY 0? IHB BSIIMAIBI) MAI'JgRIAli I B IHB ISOMIBD 
G OLIJMB ■I00II 5&; 

■/olume of concrete in the column, YOOil is given by 

VG0..1I = b^x - P) (3.1) 

Volume of steel in the column, Y.hST^ is given by the following 
expression 

VAST1 = p.b^ (Pf-D) (3.2) 

Volume of steel used in the footing can be estimated by 
using the following expression, 

■ VA3T2 = .AST2 x I. t .1313 x B (3.3) 

Volume of concrete in the footing has two components. Pig, 2. 2. 
(a) Volume of portion A and (b) volume of portion 3, 

Volume of the rectangular portion . ( A) is given by. 
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VI = 0.15 X Lx B 


and volume of remaining portion (3) is given by, 


V2 


1)1 B-b 

/f B 


,L-d 

1 - (“^)y)]dy 


On integration, 
31 

72 = 


(2.L. 3 + B.b + L.b + 2b^) 


(3.4) 


(3.5) 


where, D1 - i)--0,l5 m and 0.1 5 m being the value of the 
thickness of the footing at the edge. 

.'low, 7001^2 = V1 +• V2 ( 3 . 6 ) 


where , VOO:T2 is :he volume of concrete in the footing. 

Ihus total volume of concrete and volume of steel in the 
entire isolated column footing are given by 

VCOV = ¥00^11 + 700112 (3.7) 

and VAST = VA3T1 + VAST2 (3.8) 

respectively. 

Volume of excavation (VOX) and volume of back filling (VTIi) 
are estimated as follows: 

VOX = Df.L. 3 . (3.5) 

and Vj?IL = VOX - VCOX (3.10) 

5*3 GOST 0? MATER IALS: 

The prevailing market rates are shown in Table 3*1» 
These rates have been used in the analysis to estimate the 


total cost. 
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3*4 SMMj ME 0? IHji! PR03mi; 

J'ig, (2.2) shows the plan and sectional elevation 
of an isolated column footing. G-iven the loads and moments 
on the column and the data regarding the soil profile, the 
problem is to determine the design variables, i.e., the 
dimensions of the column, dime.nsions of the footing and the 
depth of embedment in such a way that the total cost of 
the isolated column footing is a minimum. 

Analvsis is to be carried out when the soil test 
results are available from common type of field tests (e.g. 
static cone penetration test, stemdard penetration test) and 
standard laboratorjr test (e.g. standard tri axial test and 
oedometer test). 

3.5 MAL YSIS : 

3 . 5 • 1 Des ign Variables: 

Die following are the six design variables which 
control the cost of the isolated column footing: 

(1) Breadth (h) of the column, 

(2) Percentage of steel (p) in the column 

(3) Depvh of embedment -(Df) 

(4) .Breadth of the footing 

(5) Ratio of length to breadth of the footing (m),and 

(6) Diickness of the footing (D) at the centre. 
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5*5*3 ObJ e G ti ve Pune ti on ; 

Ihe total cost of the single isolated column footing, 
which includes the cost of concrete, cost of steel, cost of 
ejccavation and cost oi hack-filling is taken as the objective 
function. 

Ihe objective function '5'' is given by 


? = C1 X VCOil + X:2-X TASI + 05 x T'S + 04- x ‘Vl’IIi - 

(5. 1 1 ) 


where, C1 = cost of concrete 

02 = cost of steel 

03 = cost of excavation I 

I 

. i 

and 04 = cost of back filling. I 

7001 = volume of concrete in the isolated column footing | 
7ASI = total volume of steel in the isolated coluimn footing j 
70X and 7FIL are the volumes of excavations and volume 
of backfilling respectively. ! 


3.6 DSS I (11 005S TRAI XIS : 


Ihe following restrictions are imposed on the design 
variables so as to ensure that the optimum design is a feasible 
design from the engineering point of view . 

1 ) UX X 3/6 (3*12) 

2) H Y < 1/6 (3*13) 

where, ;SX and BY are the eccentricities in the x and y 
directions respectively. Ihis .ensures that the resultant lies 
with in the middle third of 

5 ) ^ax '^all (3.14) 


sue ^ umsk 

. ^ . 

, A 17180 ■: 
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Where, ^nax maximum soil pressure a-r y^ is estimated by 

using equation (2.10) and qg^]_ is the allowable bearing 
pressure of the soil and is estimated by using the equation 
(2.44) . 

4) 0.0080 < p < 0.06 (3.15) 

'Ihe I.S. 456-1978 specifies that the percentage of steel (p) 
in the column should lie between 0.80 and 6,0 percent. 

MUX a MUI a 

5) ( ) + < -i^o (3.16) 

MUX1 MUY1 


6) > W 

where P^ is the ultimate load of the column given by 
3reslar, (Ref. Mallick and G-upta, (1979)) as 


u 


P P . P 

f_ux * -^uy ' ^ 0 _ ___ 

^ux * ^0 ^uy‘^o“^ux * ^uy 


(3.17) 


ux 


uy 


ultimate load under eccentricity, EX (SY = 0) 
ultimate load under eccentricity, EY, (BX =0) 

Pq = Ultimate axial load capacity. 

7) Dfi < :0f < 3 (3,18) 

wnere, = mi?iimum depth of embedment (meters) required given 

by the equation, (jain and Jai Krishna, (1973))j 


1)4 


1 


P 1 - Sin 0 p 

^ 

^ 1 + Sin 0 


(5.19) 


2 

where , P — pressure in kg/m on the soil 

3 

) = unit weight of soil in kg/m 

0 = angle of repose of the soil. 

8 ) < 4 ^ ' 


( 3 . 20 ) 



where, transverse shear given by equation (2.24) and 

Ca ~ a<llowable shear stress in Ml 5 concrete given by 
equation (2.25) 

9) <lv < (3.21) 

where q^ = shear stress at a distance ’d’ from the face 
of the column given by equation (2.28) and q„ = allowable 
shear stress conforming to IS : 456-1978, which depends upon 

size of the foundation and area of longitudinal reinforcement. 

10) ?or safe loisd transfer, the bearing stress 
should satisfy the following condition, 

tW 1 0.45 fck fWAi (3.22) 

1 1 ) Jbtal settlement ’Sg ' should be less than or 

/ 

equal to the permissible value. 

Thus for footings on clays, Sg £ 65 mm (5.25) 

and for footings on sands, Sg < 40 mm (3.24) 

1 3 ) Differential settlonent should be less than 
equal to the permissible value 

i.e. DI?? < 0.0055 -x Span (3.25) 

where Dlil' is the differential settlement between any two 
points and Span is the distance between those two points. 

14 ) :?actor of safety (S-?) agcainst overturning 
moment should be greater than or equal to the specified value. 

i.e. S. 7. >1.5 (3.26) 

and , 

15 ) lector of safety against sliding should be greater 
than or equal to the specified value. 
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The notations used in this chapuer are the same .as 
those used in the previous chapters. The total number of 
constraints in the present problem is 25. 

3.6. 1 R esults and Sla cuseions : 

A software program has been developed for the settlement 
controlled opt-irnal design of footings. All computations were 
carried out on Dec 10 computer system. To study the effective- 
ness and correctness of the computer program, limited results 
on the various aspects of the present study have been obtained 
and are presented as follows : 

it is a well known fact that efficiency of any numerical 
scheme depends to a great extent on the initial starting point. 
As such, to study the effect of the starting point on the 
final solution, various starting points were • taken. Such 
studies are also essential to know whether the obtained 
solution is a global optimum or not. 

'Table 3.2 gives the optimum value of the objective 
function obtained when tried with different starting points. 

It can be concluded from the table that the starting point 
does not have significant influence either on the value of 
the design variables or on the value of the objective function. 
However, it has been found that selecting the starting point 
closer to the optimum leads to lesser number of function 
evaluations. Yariation in the p.P.H. time for different 
starting points is very small. 
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One of the shortooraings of the penalty method is the 
severe illconditioning of the unconstrained minimization 
problem as the value of the penalty parameter tends to its 
liraiting value (infinity or zero). To study this aspect and 
to see whether minimization of 0 leads to the minimization 
of jp , results were obt,ained and presented in j?ig. 3.1. The 
figure demonstrates that as the solution procedure progressed, 
when the number of function evaluations is aoout 1800, the 
minimization of 0 leads to the minimization of 1. Hence 
it can be concluded that there has not been any ill conditioning 
and the scheme works well for this class of problems. 

Studies were carried out to check the convergence of the 
sequential unconstrained minimization technique when applied 
to this problem and the results are presented in T’ig. (3*2). 

It can be observed from the '.lig. (3.2) that as uhe penalty 
parameter (]?) decreases •'’rom ^0 to 0,1, there is a sharp 
decrease in the objective function value. As 'r' is further 
decreased, the decrease in objective function is not apprecia- 
ble and the solution converges to the optimum value very slowly. 
Stresses below the foundation have been calculated by dividing 
the foundation in to 100 strips. If the dimension of each 
■strip is greater than 1/3 of the depth at which the stress is 
calculated, then the stress at that depth is calculated by 
means of linear interpolation of the stresses between the 
surface gcd the next depth. However, no interpolation is 
done if the depth at vtoieh th©. stress is required is greater 
than 3 times the size of e^h effect of the aoove 
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mentioned apnroximation on the final solution has been studied 
by dividing the foundation in to 225 strips, fable 3.3 gives 
the design variables and friction values for both 100 and 225 
strips. ]?rom tne table 3.3) ii can very well be concluded 
that final solution is not at all significantly affected 
by the approximation. However, the O.P.TJ. time required 
for 225 strips is about 3 times more than that for 100 strips 
without yielding any significantly better solution. 

Studies were also undertaken to check how much net 
saving one may make with the help of the computer aided optimal 
design as enunciated. After doing several trials, an economic 
and one of the best possible designs which may be obtained 
mfanually is done and the dimensions so obtained by this design 
is selected as the starting point for the optimization 
procedure. Ihe optimum values have a sig’iif leant saving over 
this initial design. To illustrate this, the following two 
problems were solved by the. developed program and the results 
are summarised, 

3.7 I LLUS'TRAIIOi; 

Probl^ _ 1 : 

load on the column = W = 750 kil 
Moment about X-a.xis, MUX = 100 kl-m 
Moment about Y-axis, MUY = 100 kX-m 
Horizontal load about X axis. Hi = 6 kX 
■ Horizontal load about Y axis, H2 = 6 kX 

M1 5 coTicrctB 33 id ^250 s'fcoGl iiavo bee/i usod* 
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Parameters of the soil are, 
average G = 27.0 kl/m^ 
average 0 = 24?6 
average corrected IT value =10 
depth of water table below the footing = 1.0 m 

Bie fable 3. 4 gives the cone-penetration v.alues at 
different depths, 

'Ihe results and the percentage savings over the 
initial design are presented in - fable 3.5. 

Probl em 2: 

Problem 1 is analysed as a layered soil. 

Coefficient of volume change, *mv' for each layer is computed 

from the relation, 

1 

mv = — — 

f2 X H 

2 

wnere , mv = coefficient of volume ch^mge in m /k]7. 

Because of the non-availability of data regarding 
undrained Young’s modulus (E^) for each layer, an average 
value of equal to 10000 kl/m , has been assumed for the 

present problem. Talue of f 2 been taken equal to 700, 
Equations (2.6S) and (2.5?) have been employed to calculate 
the initial and consolidation settlements respectively. 

'Ihe results and percentage savings from the initial 

design are given in Table 3.6. 
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5* 8 G!-:5im il, DIS (yjss 10 : 

Hie general observaliions made from "the observed 
resulljs and i bs probable regi ons are discussed and summarized 
as follows: 

i) Ihe optimum is reached when the equation (2.6) is 
almost crltioally satisfied. 

ii ) Ihe depth of embedment 'D^* has direct or indirect 
influence on the value of the objective function in the 
following m.mner: 

a) Increase in the depth of embedment increases the 
allowing bearing preseure significantly. Higher the value of 
the depth of embedment (D^), the lower is the plan area of 
the footing required. 'Ihis in turn leads to a considerable 
saving in the area of concrete req>uired in the footing. 

b) Dotal settlement also decreases with the increase 
in the depth of embedment (Sf)- Dhe settlement, and not the 
bearing capacity is a major factor which controls the required 
plan dimensions of the footing. Besides, the decrease in the 
total 3 ettl( 3 ment also decreases the differential settlement 
between the different points in the footing. Dhus the decrease 
in settlement contributes to the saving in, the plan area of 
concrete in the footing. 

c ) 5*0 omi jigs are never placed on the ground sun ace . 

They are always placed at a depth i'D^) below the ground surf ace. 
?ox (J943), has given the coefficient (D'3) wnich should be 
multiplied with the computed surface settlement so as to obtain 
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the actual settlement at that depth, Ihis coefficient is 
also lound to decrease the settlement with the increase in 
iD^/ 3 ) ratio, fliis in turn reduces the actual settlement 
which has an effect over the plan area of the footing as 
discussed in (b). 

All the factors (a), (b) and (c) contribute towards 
the saving in the plan area of the footing, 

iiowever, the increase in the depth of embedment (D^) 
increases the volume of concrete and volume of steel in the 
column. Besides, the moment coming on the column due to a 
horizontal lo.^ on the column also increases with the increase 
in the depth of embedment (B^). But the saving in the volume 
of concrete in the footing is relatively much higher than the 
extra quantity of concrete used in the column. It should also 
be noted that the depth of embedment (3^) increases the 
allowable bearing pressure to a maximum value equal to, twice 
the allowable bearing press'ure at the surface. Ihe extra cost 
involved due to the increase in the depth of embedment (L^) 
on the cost of excavation and cost of filling ;are relatively 
low when compared with the cost of concrete saved in the footing. 
Consequently, 'che present work enables to select the optimum 
value of depth of foundation (Df) required for a pa-.-ticular 
soil, and for a particular foundation load so as to save the 
area of concrete in the footing as much as possible. 

iii) It has been observed that differential settlement 
between any two points in a footing actually, controls the 
minimum dimensions of the fO;Oting,'^e4iiired so as to keep the 
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ditferential settlement with in the allowable limits. This 
aargument is better .justified especially when a foundation is 
subjected to eccentric locjds. 

iv) hie overall depth of the foundation is adjusted 
along with the plan dimensions of the footing in such a way 
that volume of concrete used in the footing is a minimum, 
Bvisides, area of steel in the footing will also be taken in 
to consideration while adjusting the value of '0. 

3-9 QQ;- CTI3I Q^Si 

( 1 ) \ so.f tware for the op Cimal limit state design of 

an isolated column footing subjected to bi-axial moment 
has been developed. t 

(2) lOie program developed is a general one and can handle 
both homogeneous and non-homogeneous soil profiles. 

(3) Input for the soil-data can be the test results of 
standard penetration test, cone penetration test and 
triajxial tests. 

(4) Hie user is only required to feed the soil data, forces 
on the column and the characteristic strengths of steel 
and concrete used, 

(5) hae developed program ensures atleast 8-10 percent 
saving in the totiaL. cost of the structure. 

GJhe program will be available for reference. 
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(EABIE 3.2 I EEFSGT 0? STAR'EnTG POINT ON THE OBJEGTIYB PGNGTION 
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table 5.3 ; 


EEEBGT 0? SIZE OE STRIP OIT THE OPTIMUM SOLUTIOH 


SI. No. Design variables Number of 

- — strips 


1 0.450 0.0142 2.01 3.45 1.008 0.635 100 


2 


0.452 


0.0137 2.00 3.44 1.010 0.631 225 




TASItS 3.4 s JEST R^SSULK OF OOSIE PEHSTRATIOIT TEST 


Depth, in. metres 

Gone penetration 
value in 

Thickness of 
each Iyer in 
metres 

0.5 

2500 

1.0 

1.5 

2800 

1.0 

2.5 

3200 

1.0 

3.5 

5600 

r.o 

4.5 

3900 

1.0 

6.0 

4200 

1.0 



I ‘3 Li. '5.5 


CP rii LJr( J3S la:? 7 IR: 1 3 J 5S OJ P:‘; IOR JLIG,,:; 
S 171.71 15 COST- :.I0;. 105.3 :IJ0 US SOIL 


51, 7alue nf the '/^iriablGs 

• 

b p .3 m-'fy'3 


?unction 

^Tialae 

(cos'o) 


o / 0 

.Savings 
i'l Jost 


1, Initial 0.50 0.0100 

2 . Op ';irrxi 0 , 4.5 0 . 0 1 12 


■ts* ■' , . - 

1,S0 3.65 0.955 O.^SO- 5076.97 

10.93 

2-01 S.45 1.008 0.63-^ 4576-i1 
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Function values (F and ^ 


D(1)= 0*55 
D{2)= 0*096 
0(3) = 1*2 
D(4)= 6*0 
0(5) = 1*206 
D(6)= 0*50 



Number of function evaluation 


FIG-3-1 PATH OF F 

POINTS 


and <P functions FROl 



value (F) in thousands 


STARTING POINT 


1Bh 


D{1)= 0'55 
D(2)- 0'096 
0 ( 3 )= 1*2 
D(4)- 6-0 
D(5)= 1*206 
D(6)= 0*50 


«|5-FUNCT10N 





F-FUNCTION 


10 ^ 10 
Penalty-parameter (r) 


ATION OF F WITH PENALTY PARAMET 


GZAPm - 17 


OPIII"rM PLAZ DB0JZSIOi'iriZGr OP A SROIIP OP POOTIZGS 

4. ^ GSJgiAL PR IZ aiPBiiS OP BiilSIG Z: 

'Uie average ultimate bearing capacity of the footings 
in a group is greater than the ultimate bearing capacity of an 
isolated footing of the same size. ^Wien the footings are 
placed in a group, the individual foundation failure mechanism 
becomes distorted due to the overlapping of the failure surfaces. 

'fhe settlement for a given load intensity decreases as 
the centre to centre spacing between footings decreases, Singh 
et al. ( 1973 ) . As the footings sink in to the ground, the soil 
compresses to a certain distance on both sides of the footing. 
This results in an increase in the strength of soil. 'The 
increase of strength of soil plays its role in the settlement 
characteristics, Singn, et.al. (1973) . 

'The ultimate bearing capacity and settlement of ad;iacent 
square footings on sand are not significantly affected when 
the spacing between them is more than five times the width of 
the footing, Singh et.al (1973)- 

■Based on the experimental results, the following factors 
have been defined in order to account for the change in behaviour 
of a footing when placed in a group (Singh et al. ,1973). 
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(a) Interference efficiency factor for bearing capacity: 

It is the ratio of the ultimate bejaring capacity of the 

footing group to that of an equal number of identical isolated 

footings (Stuart, 1962), i.e, 
q^ (group) 

^ n X q^ (isolated) I 

where, q^ = ultimate bearing capacity of the soil, 
n = total number of fooTiings in uhe group. 

The ultimate bearing capacity of a smooth rectangular i 
footing can be predicted by introducing the interference I 

factor (l3) for bearing capacity in the Terzaghis ultimate \ 

[ 

bearing capacity relationship as follows (Singh et al.l975). I 

s i 

Pg = 2.25 - 0.31 (— ) for S/j 3 <3.25 (4.2a) i 

3 

fg = 1.69 - 0 . 135 (- 5 -) for S/3 >3.25 < 5.0 (4.2b) 
and = 1.04, for S/S = 5 (4.2c) i 

" I 

(b) Intei^f erence efficiency factor (ig) for settlement: I 

' 

It is the ratio of the settlement of the footing 

group at a given intensity of pressure to that of an identical 

isolated footing at the same intensity of pressure multiplied 

by the number of footings in the group. 

■S' - (.group) ( 4 . 3 ) 

n X settlement (isolated) 

■ ■ ! 

■ 

It has been observed that the factor I'd increases 

■ : 

almost linearly with the increase in S/B ratio, Singh et al.( 1973) 

' 

'Bg' is given by 
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?g = 0.4 + 0. 1 X S/3 , for S/:3 < 5 ( 4 . 4 ) 

In the abo'/e expressions, 

B = width of footings, and 

S = centre to centre distance between the adjacent 
footings, 

All the investigations on the interference between the 
adjacent footings have been done only when the footings lie 
along a line* However, no study has been made on the behaviour 
of a footing when it is surrounded by a group of other footings* 
Bierefore, in the present work, it is assumed that a footir^g is 
affected by only one of the adj-acent footings whichever gives 
the maximum 3/ 3 ratio. 

It should be noted tnat in the case of clayey soils, 
it is rather more reasonable to assimae that the properties 
of the soils are not much affected due to the presence of 
the adjacent footings. 

It is well known that the settlement generally 
governs the design of foundations on sands. Hence ignoring 
the effect of interference on the bearing capacity does not 
make the design conservative. 

Ihe present work comprises of finding out the optimum 
plan dimensions of a group of axially loaded footings for a 
given depth of foundation. 

Debeer's method corrected by Meyerhof (1965) has ‘been 
employed to find out the settlement below the footing. 
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Bie stresses and displacements caused due to load, on 
one footing below the other footing plays a significant role 
on tne shape and size of the other footings. Ihe serious 
limitation of the Schmertmann 's method of determining, the 
settlement lies in the fact that Schmertmann 's method can 
determine the settlement due to the footing load below itself 
only. Hence Debeer's method is used for calculating the 
settlement. 

4- #2 yUHQ TIOH, D ESIdH .LW OOHSmM-ifS : 

It is well known that cost of a piece of land is rising 
at a very fast rate especially in urban areas. The idea of an 
engineer should be to optimise the cost of a structure in all 
possible ways . Saving in the land area can be achieved to 
some extent by optimizing the pl.m area required for a group 
of footings. This will in addition, ensure an equivalent 
saving in the plan area of concrete. 

A group of foundations are generally designed for an 
equal settlement or an equal allowable bearing pressure . 
Gonsiderablo saving can be made by designing a foundation 
group for an allowable differential settlement instead of , 
desi.gning the foundations for an equal settlement. It is 
manually difficult to design a group of foundations such that 
differential settlement between any two foundations in a group 
is nearly equal to the maximum allowable dif-^erential settlement 
•This can be done easily in an dseful way by using the 
mathematical programming techniques. Hiis is the basic idea 
T-iAhima taking this problan as an optimization problem. 
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4.2.1 S tatement of the P:?oblemi 

G-iven the loads on each foundation, the soil 
properties and the allowable total and differential settlement, 
the problem is to determine the optimum size of each foundation 
in a group.. 

Ob.iecti've function: 

The total plan area of 'all the foundations has been 
taken as the objective function. The objective function is 
given by 

n 

f = E I. 3. (4.5) 

i=1 

where, 

= Length of ith foundation 
3i = Breadth of ith foundation 
n = i7umber of f oundationsin the group. 


4.2.2 Gonstr;a ints : 

Some restrictions are required to be placed so as to 
ensure that the foundation ;trea provided is adequate from 
the safety point of view . 


where , 


where , 


’ ') %iax^ ■- 


q 


alii 


(4.6) 


ra<aximum pressure on the soil given by 



(4.7) 


% 


a 


allj^ 


is the lotad on the ith footing 

allowable bearing pressure on the ith 
foundation given by equation (2 . 44 ) . 
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To check whether there is any illconditioning 
of the uncons tr.ained minimization problem with the decreasing 
v-alue of the penalty parameter, results were obtained and 
are presented in .?ig, 4,1. Hie figure shows that minimization 
of the 0 function leads Co the minimization of T* function and 
as such it can be concluded that no ill conditioning has 
occurred leading to any premature termination of the algorithm. 
When the number of function evaluations is about 3800, the 
minimization of 0 leads to the minimization of J function, 

■s'ig. 4.2. 

Optimum plan dime.nsioning of a group of footings is 
generally done either on the basis of equal bearing pressure 
ob? on equal settlement. Hiis approach is too conservative. 
'Considerable saving can be achieved if all the footings are 
designed in such a way that the differential settlement 
between any two footings is almost equal to the maximum 
allowable differential settlement. This design becomes 
tedious and cumbersome when the number of foomiiigs in the 
group is more. However, this kind of a problem can be 
easily hiandled with the aid of a computer, using mathematical, 
pipfgramming techniques. 

To illustrate tne amount of saving that one can achieve 
by the developed methodolog3r in present study over conventional 
approach, the following problems have been solved. 

Ill ustra,tio.n.v 

Referring to lig. (4.5) , the loads on the footings 
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71 = 800 ; '73 = 800 ; W5 = 1500 ; ’JS/7 = 800 azid W9= 800 
W2 = 1200; W4 = 1200; W6 = 1200; V© = 1200. 

Properties of the soil: 

Average G =0.0 

/ 

iverage 0 = 52° 

Averaije corrected J = 22 

Equal spacing between the footings 7.5 m. 

’Ihe cone penetration values are tabulated in the Table 4.2. 

Following the procedure as suggested by Teng (1963) , an 

initi;al design is obtained such that all the footings are 

subjected to equal bearing pressure. This ensures t-hat there 

is no differencial settlement between any two footings. The 

dimensions obtciined from this initial design are taken as the 

starting point for the optimization problem. The saving 

obtained by the optimum design over that of the initial is 

presented in Table 4.3 . It can be seen that rhere is a 

substantial net saving in 'che pliui area (41.3 percent). 

In case of clayey soils, the most reasonable asaumptiui 
we can make is that the soil properties are not changed due to 
the interference between the adjacent footings. It is known 
that the properties of cl.rrev soils are not changed like that 
of Sfand on densification. To enunciate the s.avings that one 
can make while designing a group of shal''ow footings on clay, 
the following problem has been solved by the developed 
program. 

s 

Problon 2: 

Referring to Pig. (4.3) ? the loads on each footing in kil 


are as follows : 
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Loads in lei 

il = 750 ; v/3 = 750 ; W5 = 1000 ; WJ = 750 and W9= 750 
W2 = 800 ; W 4 = 300 ; W 6 = 800 ; m = 800 
Properties of the soil are as follows: 

Average G - 21 ki/m^ 

0 = 0.0 

■^dual spacing between the footings 7.5 m. 

Gone penetration values are the same as those given in 

fable 3.5 . 'fhe average corrected, i = 10. fhe depth of 

foundation = 1,8 m. 

'fhe results and the net savings in the plan area over 
the initial mannual design are snown in fable 4.4 . In this 
case also the net saving in the total plaii area is substantia] 
( 37,3 percent). 

4 . 3.2 Influen ce of Se-pth of Poundation on the fotal Plan .ki^ei 

It is a well known fact tnat depth of embedraent (P^) 
increases the allowable bearing pressure and decreases the tci 
settlement. Depth of embedment (D^) also has a sig:nff leant 
influence on the required plan dimensions of footings. fo 
illustrate this, optimum total plan area of a ;group of footing 
are determined for different depth of embedments and the 
variation of optimum total plan area with Dp is shown in 

Pig. ( 4 . 4 ). 

It can be observed that as the depth of foundation (Dp 
increased from 1 . 2 m, to 2.0 m the net decrease in the total 
plan area (P) is quite substantial (5.5 m^) ; however , with th 



inc'ease in 1.5 m, The rate of decrease in ? decreases, 
beyond 2 m the effect of depth of foundation on the total pl<an 
area is not appreciable. The depth of foundation could have 
been incorporated in the formulation as a design Vrjriable but 
was not so considered as from various considerations , liie 
normal practice is to get the plan area for a given depth of 
foundation. However, as tne scudv indicates, this should 
also be considered as a design variable for a realistic 
analysis to arrive at the 'global optimum value as the 
solu'tion that one gets 7 !or a given depth of foundation can 
at best is a local optimum. Settlement below each footing 
is presented in Table 4»5- 

4-. 4 WJlOLUSIOHSj. 

1 . A general optimization scheme has been outlined for 

the optimal pl;3n dimensioning of a group of footings and 
it has been demonstrated that the developed algorithm is 
quite' efficient is arriving at an optimal solution. 

2. Tory significant saving in the total plan area of 

a group of footings can bo achieved by following the 
scheme. 

‘jChe obtained solution is independent of the starting 
point. 
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TABIiB 4.5 ; SBirMENI? KSICW BACH FOOTIIG 01 SAND 


Number of the 

Settlement below 

Settlement below 

footing 

the initial design 

the final design 


(mm) 

(mm) 

1 

30.870 

37.23 

2' 

51.309 

35.36 

3 

31 ♦ 020 

37.10 

4- 

31 t 324 

35.60 

^ , 

31.520 

34.40 

Ite' - 

31.328 

35.63 


31.078 

37.208 




CD 

31.297 

35.330 

9 

30.790 

37.287 






Penalty- parameter (r) 

FIG.4-1 DIAGRAM SHOWING VARIATION OF F AND 




Function value F and 
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FIG- 4-3 PLAN SHOWING THE GROUP 
OF FOOTINGS 




Tolal atQQ CF) In 
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